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Abstract 

We consider the screening of a magnetic impurity in a dr^2_y2 wave super- 
conductor. The properties of the dr^2_y2 state lead to an unusual behavior 
in the impurity magnetic susceptibility, the impurity specific heat and in the 
quasiparticle phase shift which can be used to diagnose the nature of the con- 
densed state. We construct an effective theory for this problem and show that 
it is equivalent to a multichannel (one per node) non-marginal Kondo problem 
with linear density of states and coupling constant J. There is a quantum 
phase transition from an unscreened impurity state to an overscreened Kondo 
state at a critical value Jc which varies with Aq, the superconducting gap away 
from the nodes. In the overscreened phase, the impurity Fermi level e/ and the 
amplitude A of the ground state singlet vanish at like Aq exp(— const. /A) 
and J — Jc respectively. We derive the scaling laws for the susceptibility and 
specific heat in the overscreened phase at low fields and temperatures. 

PACS numbers: 71.27.+a,75.20.Hr 
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I. INTRODUCTION 



The problem of magnetic impurities in d-wave superconductors has been a subject of 
intense research [jlj-Q. It is known that in anisotropic superconductors, such as a dx2_y2 
state, magnetic impurities act as pair breaking centers and hence reduce the amphtude 
of the condensate. Experimentally, the main effect of these impurities is to reduce the critical 
temperature Tc of the superconducting state . 

Quantum mechanical fluctuations of magnetic impurities give rise other important effects 
particularly when coupled to the fermionic quasiparticles of the superconducting state. In 
a normal metal these correlations lead to the Kondo screening of the impurity and to the 
generation of dynamical energy scales such as the Kondo temperature. In high temperature 
superconductors, the effects of the magnetic impurities appear to depend significantly not 
only on the nature of the impurity but also on where the magnetic impurity is located. The 
conventional interpretation of the role of magnetic impurities in high temperature super- 
conductors relies, for the most part, on the chemical differences of the impurities (mainly 
Zn and Ni) and on the actual location of the impurities on the lattice relative to the CuO 
planes. The main focus of recent work on this subject has focused on how much different 
impurities are able to depress T^. and on the power laws that static (or clasical) impurities 
induce on low tempreture properties. 

In this paper we investigate the physics that results from the exchange coupling between 
isolated magnetic impurities and the quasiparticles of the superconducting state. In par- 
ticular we will be interested in finding out under what circumstances there is a Kondo-like 
dynamical screening of the magnetic impurity by the quasiparticles. The mechanism that 
we have in mind is analogous to the exchange coupling between magnetic impurities and 
the electrons of a Fermi liquid that causes the Kondo effect. However, unlike the Kondo 
effect in metals, because the density of states of normal quasiparticles in a d-wave super- 
conductor vanishes at the Fermi energy, screening is absent in perturbation theory and a 
critical exchange coupling between the quasiparticles and the magnetic impurity is necessary 
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for Kondo screening to take place. Thus, magnetic impurities which couple strongly to the 
quasiparticles, such as Zn which substitute for is the planes, may actually be Kondo 
screened at very low temperatures while Ni, which appears to couple more weakly, may not 
get to be in a Kondo screening regime. Likewise, magnetic impurities on sites away from the 
CuO planes are weakly coupled and therefore are less likely to undergo Kondo screening. 
We will see below that the critical coupling (which we will only estimate very roughly in 
this work) is controlled by Aq, the size of the gap away from the nodes at zero temperature 
and tipically it is a fraction of Aq. 

The onset of Kondo screening at a critical value of the exchange coupling constant is a 
quantum critical point. We will show in this work that the behavior of the magnetic impurity, 
both near and beyond the phase transition, has unique signatures which follow from the 
nature of the condensate and hence can be used to investigate its nature. Among its most 
salient features are the temperature and magnetic field dependence of the impurity magnetic 
susceptibility and specific heat which exhibit strong deviations from Fermi liquid behavior. 
We will also show that the quasiparticle phase shift exhibits a strong frequency dependence 
and a broad resonance and that the structure of the quasiparticle scattering matrix has 
detailed information on the phases and signs associated with d-wave superconductivity. 
Thus, the physics of magnetic impurities in a d-wave superconductor can be used to diagnose 
the nature of the superconducting order. The purpose of this paper is to describe these effects 
in detail. 

The Kondo problem in metals has been intensively and extensively studied and it is by 
now very well understood. It is described in terms of a smooth crossover from a marginally 
unstable fixed point at zero exchange coupling to a stable Fermi liquid fixed point with 
a screend impurity P-[Tl|]. The validity of this picture has been confirmed by the exact 



solution by the Bethe ansatz [0,0 and by large- methods |T^. We will refer to this case 
as to the marginal Kondo problem since it bears a strong resemblance to a critical system 
at its low critical dimension. 

In a conventional s-wave superconductor the Kondo effect is suppressed by the formation 



of the superconducting gap, as shown by the classic theory of Abrikosov and Gorkov [jlS . 
However, in the case of a d-wave superconductor, there are quasiprticle states inside the 
superconducting gap which concentrate near the nodes of the order parameter. Although 
the density of states vanishes at the Fermi energy, for strong enough exchange coupling it may 
still be possible that isolated magnetic impurities may still be screend by the quasiparticles. 
The central idea of this work is that the way this screening happens may be used as a tool 
to study the superconducting order. 

Impurities cause many different effects in superconductors. In the case of a d-wave 
symmetry, any sort of scattering breaks pairs and, for instance, static magnetic impu- 
rities produce quasiparticle bound states in the gap of the superconductor. However, the 
binding energy of these states vanishes in the vicinity of the nodes of the d-wave supercon- 



ductor p,^,|T6[. Random potential scattering also leads to interesting effects, in particular 



close to the nodes of the superconductor where the density of states (DOS) of quasiparticle 
states, which behave like Dirac fermions near the nodes, vanishes linearly with the energy 
(measured from the Fermi energy). It has also been shown that random potential scattering 
should generally lead to a finite DOS at the Fermi energy (zero) for Dirac fermions in ran- 
dom potentials fT^HIHl d-wave superconductors P(]|J^,pT[] . The precise behavior of the 
DOS appears to depend on how many nodes are coupled and on what channels are mixed 



by the scattering processes P, p!9| , p2| , ^ . However, if the superconductor is sufficiently clean, 
the effective DOS induced by the disorder is exponentially small [jl^jliH and its effects can 
be neglected. Notice, however, that the combined effects of Kondo screening and random 



scattering is a problem that is still not undertood, even in metals. [24,25 



The problem of a quantum magnetic impurity coupled to fermions with a vanishing DOS 



at the Fermi energy was first discussed by WithofF and Fradkin |^6[. In contrast with the 
Fermi liquid which has a finite DOS at the Fermi energy , Kondo screening of the magnetic 
impurity could only happen for values of the exchange coupling constant J larger than 
some critical value. If the DOS vanishes with a power of the energy with exponent r, the 
fixed point at J = is stable for r > and a new unstable fixed point appears at Jc > 



signaling a quantum phase transition. We will refer to this problem as the non-marginal 
Kondo problem and r measures the deviation from marginality. Using a close analogy with 
the theory of critical phenomena, Withoff and Fradkin developed a large-A^ theory for this 
problem and found that the essential singularity of the Kondo temperature is replaced by 
power law singularities determined by the DOS exponent r. However, in that work only 
the regime < r < | was considered. The behavior of these systems for small r was also 
exlored by Chen and Jayaprakash and by Ingersent who used a generalization of 



Wilson's numerical renormalization group (RG) for this problem. The case of interest for a 
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d-wave superconductor is r = 1 which turns out to be special in several ways. In ref. [P9[ 
discussed recently the closely related problem of a magnetic impurity in a flux phase which 
is also an example of an r = 1 system. In this paper we show that these two problems can 
be mapped into each other in spite of the fact that spin fluctuations do break Cooper pairs. 

One aspect of the physics of magnetic impurities in a d-wave superconductor that we 
will not consider here are the effects of the depression (and/or actual vanishing of) of the 
d-wave order parameter near the impurity site. This effect is independent of the spin and 
it actually static. In any event, the vanishing of the condensate at the impurity site leads 
to terms that do not conserve fermion number in the effective Hamiltonian and hence may 
lead to important (spin-independent) effects. This is a conceptually important problem 
and it will be addressed elsewhere HSOll. A self-consistent calculation based on the BCS 



approximation can be found in ref. [23 



In this paper we will make use of a very simple model of the quasiparticle dynam- 
ics in a dx2_y2 superconductor ||3T|j3^ . We will use the fact, strongly supported by the 



corner-junction interference experiments as well as by angle resolved photoemission 



spectroscopy (ARPES) p^pSl that the high temperature superconductors have a dx2_„2 
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condensate with four symmetrically arranged nodes where the quasiparticle gap vanishes. 
The first evidence that the gap vanishes at the dx2_y2 nodal line was reported by Shen et.al, 
. Ding et.al, reported measurements of the momentum dependence of the supercon- 



ductor gap in Bi2Sr2CaCu208+z consistent with a gap function of the form cos(A;xO —cos{ky), 



as expected for a d-wave order parameter. Interestingly enough, in underdoped systems, 
photoemission supports the idea that the gap may survive through a large range of tem- 



perature into the normal state p^j35|] . It is also well established that the high temperature 
superconductors are not conventional BCS systems in the sense that their normal states 
deviate strongly from the predictions of Fermi liquid theory and that the interactions are 
strong. Thus, a straightforward BCS self-consistent approach should not work, particu- 
larly in view of the fact that there isn't a well established mechanism for superconductivity 
in these materials. Nevertheless, whatever actual the mechanism is, it should describe a 
system with nodes and with gapless quasiparticle branches. The actual coefficients of this 
effective hamiltonian cannot be derived in a simple minded way from a microscopic system 
but its form will be determined by the requirements of dx2_y2 symmetry. Thus we will use 
a phenomenologically-motivated BCS-like model for the quasiparticles with nodes consis- 
tent with dx2-y2 symmetry but without a self-consitent derivation of its coefficients. We 
will consider the case of a very clean system and at very low temperatures so as to neglect 
fluctuations of the amplitude of the superconducting order parameter. 

In section |T| we derive an effective Kondo-like hamiltonian for the problem of a single 
magnetic impurity in an otherwise perfect d^i-yi superconductor. In this model we focus on 
the effects of the quasiparticles close to the nodes of the dx2_y2 state within an energy range 
Ao, the gap of the superconductor away from the nodes. We also include, albeit in rather 
crude fashion, the effects of the states above the superconducting gap since they affect the 
value of the critical exchange constant. By expanding the electron operator in the exact 
quasiparticle states of the dx2_y2 superconductor, we map this problem of two-dimensional 
physics into an effective one- dimensional system of chiral fermions coupled to the impurity. 
The effective Hamiltonian is almost identical to the problem of a magnetic impurity in a 
flux phase that we discussed in reference ||29| . The only difference here is that the symmetry 



is SU (2) (spin) but there are four species (or flavors) of chiral fermions, one per node. 
The mapping of the electron operator into the effective one- dimensional fermion contains all 
the information about the coherence factors of the dx2-y2 state and, hence, it includes the 



pair breaking effects caused by tfie spin fiuctuations of the impurity. As a bonus, we get 
explicit (although qualitative) relations between the effective coupling constants, the relative 
importance of intra-node and inter-node scattering processes and important parameters such 
as the location of the impurity (relative to the CuO plane) and the superconducting gap 
Aq. We will find that, in fact, there is always one effective channel that matters. 

A considerable number of theoretical tools have been developed to study Kondo systems. 
For magnetic impurities in a metallic host, which have an essentially constant density of 
states near the Fermi energy, the different methods complement each other in a manner 
that we now have a rather complete understanding of this phenomenon at a non-perturbative 
level. However, with the exception of large-N methods or Wilson's numerical renormalization 
group, all the other methods (including the powerful mapping to one-dimensional logarithmic 
gases, the exact solution via the Bethe ansatz and the conformal field theory approach) 
cannot be applied to systems with a vanishing density of states. For these reasons in this 
work we use the large- A^c approach, even though N^ — 2 for the d-wave superconductor. In 
conventional Kondo systems Nc = 2 and Nc = oo are known to be smoothly connected and, 
although it is likely that this will also hold for a d-wave system, there is still no evidence 
that it is also true for this problem. In any case, given the lack of alternative approaches, 
we will present here a large- A^c theory of our problem. 

Our large- A^c theory predicts that magnetic impurities in clean cuprate superconductors 
should undergo a quantum phase transition at a critical exchange constant whose typical 
value is crudely estimated to be below the superconducting gap Aq (details are given in the 
next and in the last two sections). Our estimates indicate that, for an exchange coupling 
J with strength about 10% larger that the critical coupling, the anomalous behaviors that 
we predict should be accessible to measurements of the low temperature heat capacity and 
magnetic susceptibihty (such as in NQR) with magnetic fields ~ 1 — 10 Tesla and at 
temperatures T ~ 1 — lOK. The magnetic fields should be in-plane so as not to disturb the 
kinematic properties of the quasiparticles on the CuO planes. At temperatures higher that 
Tk (but still below Tc) the systems behaves as if it were at its quantum critical point at Jc- 

7 



In ref. we discussed the solution of the problem of a magnetic impurity in a flux 
phase problem which, as we indicated above, is closely related the problem of a d-wave 
superconductor. In that work we used the large-A^^c limit to investigate a similar phase 
transition. However, the results that we present here for the scaling behavior of the physical 
observables near this transition disagree with our previous work of ref. |2^. The reason for 
the discrepancy is that in the process of carrying out the large-A'^c limit divergent series 
need to be handled and in ref. these series were regulated in in a manner that is 
incompatible with an integer filling Q/ of the impurity in the slave fermion representation. 
This inconsistency is removed in this paper and the results that we present here supersede 



those of ref. [2S 



This paper is organized as follows. In section H we present a summary of the main 
results of this paper including experimentally accessible predictions for magnetic impurities 
in high temperature superconductors. In |ITl| we describe the mapping of the model of a 
single magnetic impurity coupled to a dx-^_y2 superconductor to an effective theory of chiral 
fermions in one-dimension with a non-marginal coupling. In section |rV| we discuss the large- 
Nc approximation and in section |V| we discuss the solution of the saddle point equations, 
valid in the ^ oo limit, and use it to investigate the phase diagram of this problem 
at zero temperature and zero magnetic field. In section ^ we calculate the low and zero 
temperature magnetic susceptibility of the impurity at zero and finite (but small) fields in 
the Nc ^ oo limit. Similarly, the impurity entropy and specific heat (at low temperatures 
and fields) are calculated in section [VI I| again in the oo limit. In section [VIII 



we 



conclude with a discussion of the implications of our results and their relation with other 
work, particularly the RG work of Ingersent. Relevant details of the computation of various 
integrals are given in the Appendix. 
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II. SUMMARY OF RESULTS 



In this section we give a brief summary of our main results and discuss their imphcations 
for magnetic impurities in high temperature superconductors. 

1. We introduce a model for a single magnetic impurity in a d^^^yi condensate (see section 
lIT|) with Hamiltonian H 



d^x J{x) 4(x)wc<^/(f) (2.1) 



k.cr 



This is accurate at energies and temperatures low with respect to the gap Aq of the 
dx2_y2 condensate away form the nodes. We use a simple lattice model for the d^-^^y^ 
quasiparticles (see ref. [|l|), which is only a cartoon of a realistic superconductor, 
but it has the correct nodal structure and that is all that we actually need to know. 
Next we construct an effective model of one- dimensional chiral fermions coupled to 
the impurity: 



H, 



dp 



= E E 

+ E E iJa/2) 
a=l a,v=-\,l 



— VVV' P dl^{p)daa{p) 

dp 



oo 

-oo 27r 



WLip) 



dp' 
2^ 



\p'\dau{p) 



(2.2) 



where v and v' are the velocities of the quasiparticles of the d-wave state along the 
two lattice directions. This is the model that we actually investigate. The nature 
of the dx2_y2 state is present: (a) in the momentum dependence of the interaction 



between the chiral fermions and the impurity (see Eq. p.32| ) and (b) in the number 
and angular momenta of the channels that are coupled to the impurity (see Eq. |3.23| 
and Eq. |3.24| ). We truncate the momentum dependence of the interaction beyond 
a momentum scale ~ Aq/(2tc\^vv') where is saturates indicating that the effective 
density of states is nearly constant for states above the gap Aq and up to an upper 
cuttoff scale D, the bandwidth of the normal quasiparticles e{k). A typical value of 
Aq for a CuO superconductor is lOOi^' and D/Aq> 10. 
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2. In sections |V| and we solve this effective model in the large Nc limit. The physical 



properties of the system in the N^. ^ oo limit are determined by the behavior of 
the phase shift 6{e) of Eq. [4.8| . This is the phase shift acquired by the quasiparticles 
of the superconductor as they scatter off the magnetic impurity. We find that, in 
contrast with the conventional Kondo problem, the phase shift 6{e) has a strong energy 
dependence. 6{e) is parametrized by the singlet amplitude A and the impurity Fermi 
level e/ (which plays the role of the Kondo scale). These are determined by solving 
the saddle point equations (Eq. |5.1| and Eq. |5.2|) which yield A and e/ as functions of 
the impurity filling Qf, the exchange constant J, the temperature T and the magnetic 
field H. 

3. Phase Transition: We find that , at T = H = 0, the system has a quantum phase 
transition at a critical coupling constant Jc ~ (2ff'/Ao)/(l + ln(D/Ao)). This transi- 
tion separates a weak coupling free phase in which the impurity is nearly free with a 
Curie-like susceptibility, from a strong coupling phase where the impurity is screened. 
We find that the impurity is actually overscreened since the impurity magnetic sus- 
ceptibility vanishes for J>Jc at T = if = 0. A naive use of the BCS estimates yields 
Jc ~ D sin^(A;o)/(l + ln(D/Ao)) which is larger than Aq unless ko < 1 (which is not 
unreasonable). However, strong coupling corrections in the superconductor (which, at 
best, can only be estimated and depend strongly on the details of the mechanism of 
superconductivity) in general will invalidate the simple relation between the velocities, 
fco and Aq. However, we expect that Jc < Aq, typically being of order of Ao/2 or so. 
However, a precise estimate requires a more sophisticated calculation than the one we 
do here. For the puroposes of this work it will be sufficient to know that < Aq. 

4. Kondo Scale: We find that, close to and above Jc, the Kondo scale = e/ is related 
to the singlet amplitude A by (Eq. |5.9| ) 



I 2x 1 \ 

e,{x. A) = Ao exp ( -- ^ _ ^ ^ + ^ _ ^ J [1 + 0( A, A In A)] (2.3) 
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where e = 2.7172 . . . and x = Qf/Nc- This singular relation between e/ and A results 
from a logarithmic singularity in the saddle point equation for the impurity occupancy. 
This singularity is absent for r < 1. In this sense, r = 1 is like an "upper critical 
dimension" for impurity problems. The actual dependence on the coupling constant 
is determined by solving the equation of state Eq. |5.15| 

1 1 / Ar 



-^)a + 0{A') (2.4) 
nvv' / 



Jc J \ nvv 

A rough estimate of Tr- can be obtained by setting Qf = 1 and = 2. We find that, 
for (J/ J,) - 1 ^ 0.2 and Aq ^ lOOK, Tk = ej ^ 9K, while , for (J/JJ - 1 ^ 0.1, 
Tk ^ IK. 

5. We have calculated the impurity susceptibility and specific heat in the overscreened 
phase J > Jc for magnetic fields and temperatures H,T < Tk which can be realized 
unless J is too close to Jc- For T <^ H <^ Tk we find (Eq. |6.15|) the susceptibility 



Ximp ~ Nc if In^ while , in the opposite regime H <^ T <^ Tk (Eq. |6.16|) 



.2 

Ximp ~ 2A^cln2 (^-^^ Tln(^). Similarly, the specific heat in the regime H <^T <^ 
Tk is (Eq. p|) Cimp(0,T) ^ %{'i)Nc^T'^\a{^), while for T < if < we find 
instead the result (Eq. ^ Ci^p{H,T) ^ N^^ (tf H\n{j^). The low field 
regime is clearly very different from a Fermi liquid although a Wilson Ratio can still 
be defined and it is finite (Eq. |7I6| ) ^'"'"'■^^^^ !^ The behavior in the high field 

regime is more like a Fermi liquid. These behaviors should be accessible to experiments 
in clean samples of cuprate superconductors at magnetic fields of 1 — 10 Tesla. 

We have not investigated yet the quantum critical regime H,T > Tk, which we will 
discuss in a separate publication [RO]. 



III. THE MODEL 



In this section we construct the model that describes the coupling of the quasi-particles 
of a d-wave superconductor to a localized magnetic impurity. We will show explicitly that 
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this model maps exactly onto a model of a magnetic impurity coupled to the spinous of a 
flux phase that we discussed in ref. I^Qj. The strategy that we will follow in this section 
consists of first writing down a simple model for the quasi-particles of the d-wave super- 
conductor with a physically reasonable coupling to a local magnetic moment. Next we will 
carry a dimensional reduction of this problem down to a model of effective one-dimensional 
(chiral) fermions which has all the symmetries of the d-wave superconductor. The effec- 
tive one-dimensional model coincides exactly with the non-marginal Kondo problem that 
was discussed in reference [^. In the next section we will use the results of to draw 



conclusions on the effects of magnetic impurities in d-wave superconductors. 

A. Free Hamiltonian 

We begin by choosing a model of a d-wave superconductor with the form of a BCS-type 
Hamiltonian. It has a kinetic energy term (which we choose to be of the form of a tight- 
binding Hamiltonian) and a pairing term with d-wave symmetry. In subsection |11I B 



we 



will describe the way magnetic impurities couple to the quasi-particles. Here we will make 
the phenomenological assumption that there is d-wave pairing regardless of the mechanism 
that gives rise to that pairing. BCS-type models which exhibit d-wave pairing (driven by 
antiferromagnetic fluctuations) have been proposed by Bickers, Scalapino and White ||3^ 



and by Monthoux and Pines [l]]. Here we will use a BCS model of this type to describe the 
dynamics of the quasi-particles. 

What will be important for the dynamics is that the model exhibits four nodes where the 
gap vanishes and that the gap is fairly large away from the nodes. Thus, we will concentrate 
on the behavior of the quasi-particles close to the nodes. Instead of using the full detailed 
form of the gap, we will replace it by a linearized spectrum with a wavevector cutoff A (rela- 
tive to the location of the node) such that the energy of the quasi-particles with wavevector 
A is approximately equal to the value Aq of the superconductor gap away from the nodes. 
The actual structure of the quasi-particle spectrum away from the nodes will play very little 
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role and such states will be neglected. This view is supported by recent photoemmission 
experiments by Shen et.ai, ||3^ and Ding et.ai, in YBaCuO superconductors where 
a Fermi surface is seen at optimal doping and it disappears progressively away from the 
nodes (where the gap vanishes) for underdoped systems. Hence, the important features of 
the quasi-particle spectrum that we will keep are the four nodes where the excitations are 
gapless, the correct behavior under lattice symmetries (and parity) and the Fermi velocities 
at the nodes. 

The Hamiltonian for the quasi-particles of a BCS-type superconductor in the absence of 
impurities is 

^0 = E e(^)4c,^. - E mc\fUi + (3-1) 

k^(7 k 

To make the model concrete we use a lattice model for the quasi-particles with a bare energy 
e{k) of the form 

e{k) = t{-k) = -2t (cos(A;i) + cos(A;2)) + At (3.2) 

The Fermi operators ct create quasi-particles with momentum k, spin a and energy e(^). 
Here fi is the chemical potential for the quasi-particles. The gap function A{k) for a d^2_y2 
superconductor given by 

A(fc) = Ao (cos(A;i) - cos(A;2)) (3.3) 

Here Aq set the scale for the gap away from the nodes. The one-particle spectrum of this 
simple Hamiltonian agrees qualitatively with the observed photoemission spectrum. 
Next we write the quasi-particle operators in the Nambu-Gorkov form |3^ 



m = I 1 (3.4) 

In terms of the Nambu-Gorkov spinors, the free part of the Hamiltonian Hq can now be 
written as 

Ho = J2 "^^k) [[e{k) - /i) ra - A{k)n] <l>{k) (3.5) 
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where ti and are two Pauli matrices. This model has four nodes at the points in the 



Brillouin zone given by (±A;o,±A;o), with ko = arccos(yu/4t). The spectrum of the quasi- 
particles crosses the Fermi surface at those points and the gap closes. As a consequence, the 
quasi-particles have a linear dispersion relation in the vicinity of these nodes. This can be 
shown 10 by expanding for small momentum departures around the (i/co, ^k^) points. 

The next step in the construction of an effective low-energy model is to describe the 
dynamics of the quasi-particles close to the nodes. To this end we assign a label for each one 
of the four nodes. Let a = 1, 2, 3, 4 be this label and we assign the label a = 1 to the node 
{ko, ko), a = 2 to the node (— /co, — ^o), a = 3 to the node {—ko, k^) and a = 4 to the node 
(fco, ~^o)- Let (f be the momentum relative to the node. It is useful to work in the rotated 
basis pi = (l/\/2)(gi + ^2) and p2 = {l/V2){qi — ^2); with velocities v = 2\/2t sm(ko) , 
v' = -\/2Ao sin(/co), where Aq is the size of the superconductor gap at its maximum value. 
Let ^a{p) denote the (Nambu-Gorkov spinor) operator which creates a quasi-particle with 
(rotated) momentum p relative to the wavevector of node a. The free Hamiltonian now takes 
the form 



In the long-wavelength limit the Hamiltonian splits into four (anisotropic) Dirac-like hamil- 
tonians. In what follows we will refer to these four sets of excitations (which represent the 
four nodes of the d-wave superconductor) as to the four flavors (or channels). 



— 




It will prove useful for our purposes to rotate $a(^) to a new field ip, 



a 



Mk) 



1 



(1 - tn)iJaik) 



(3.7) 



V2 



and to write Hq in terms of ip, 




(3.8) 
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In Eq. ( |3.8| ) the following definitions have been used: 



e+ = \J {vpiY + {v'p2y = Vvv'p+] 6+ = tan ^ {vpi/v'p2) ; 

e_ = ^J {v'piY + {vp2y = \fvv'p^\ 6'_ = tan~^ {vp2/v'pi) 

where v = 2-\/2i(:sin(A;o), "^^ = V^^q sin(A;o). 

Next we notice the fact that, as far as the kinetic energy is concerned, p+, 6*+, 
9- are just dummy variables and that the measure in the integrals is invariant under the 
change pi,p2 into It turns out, an we show this below, that the interaction term 

is also invariant under a redefinition of the integration variables. Naturally, the quasi- 
particle operators themselves are not invariant under these redefinitions of variables. Hence, 
although all explicit reference to the anisotropy can be removed from the Hamiltonian, it 
remains quite explicit in the relation between the quasi-particle (fermion) operators and the 
fields that will describe the effective Hamiltonian, i.e., in generalized coherence factors. 

Taking these observations into consideration, Hq can be put in a much simpler form 

° -^''''-'^Um^uP) (3.9) 

Jo 2ti Jo 2n y ^^pg*^ / a=i 



where Tab is the 4x4 diagonal matrix in flavor indices diag(l, —1, —1, 1). The signs in the 
matrix Tab account for the parity of each node. Here we have only kept explicitly the flavor 
(node) indices. 

We now diagonalize the kinetic energy and expand the fields in energy eigenmodes. 

m = Mp)u+id) + ^-ip)u-{d) (3.10) 

where 

1 



y±e J 



(3.11) 



are the spinors that diagonalize the d-wave BCS Hamiltonian near the nodes. The effective 
rotational invariance around each node (in terms of the redefined momenta) enables us to 
expand in angular momentum eigenmodes around each node 
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(3.12) 



This is effectively an angular momentum expansion in elliptic coordinates around each node. 
Hq is now diagonal and takes the simpler form 

'•r^ P dp °° - 



ho = y: 



VV'p E ^ab V^i+,m(bl)V^6,+,m(bl) -^i-,m(bl)V^6,-,m(b|) (3.13) 



B. Impurity Interaction 

Now we consider the interaction term for spin impurities given by 

-f^imp = ^ ' j ^(^) ci{x)T^^ic^i{x) (3.14) 

In practice we will be interested in well localized impurities. This means that J{x) is sharply 
peaked at some point xq where the impurity is located. Realistic magnetic impurities in 
YBaCuO and other High Temperature Superconductors |^ almost always involve magnetic 
atoms which either substitute a Cu atom or hybridize strongly with it. This is the case 
for Ni which, due to its hybridization with oxygen, it is believed to behave like a S* = 1/2 
impurity spin [|I|. Similarly, Zn substitutes Cu which now behaves like a missing 5=1/2 
magnetic moment and in this sense is a magnetic impurity. In all cases of Cu substitution 
we will model the impurity as a localized S* = 1/2 moment residing at a site of the square 
lattice which we will consider as the origin. Notice, however, that O can also behave like 
a magnetic impurity in the cuprates. An O magnetic impurity sits in the middle of the 
bond instead of a corner Cu site. This case leads to more complicated form of the effective 
interaction which we will not discuss in this thesis. 

The effects of magnetic impurities on Cu sites can be modeled qualitatively in terms 
of an exchange coupling constant J{x) which couples most strongly to the quasi-particles 
at X = and decays rapidly and symmetrically around a; = 0. For simplicity we will use 
a model in which J{x) is a narrow gaussian. We can see clearly from the discussion that 
led to the effective free Hamiltonian, that the only properties of J{x) that are important 
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are the amplitudes of its Fourier transform at the relative wavevector of the nodes. These 
amplitudes play the role of the effective coupling constants. Physically, the strength of the 
exchange coupling is determined by an overlap integral which decays very quickly. Thus, 
impurities which substitute Cu atoms in the plane are more strongly coupled than those that 
substitute Cu out of the plane. Also impurities on sites other than Cu sites are more weakly 
coupled to the quasi-particles than those on Cu sites. These observations are important since 
we will see in section that the impurities are Kondo screened if their exchange coupling 
constants are large enough. 

We now proceed to find the contribution of the impurity interaction to the effective 
Hamiltonian. In momentum space Eq. (|3.14|) becomes 



(27r)2y (27r) 
In terms of the NG spinors it reads 

i^imp = LjlUik- k') \S^ + S_ e,,^\{k)^]{-k') + h.c] } (3.16) 

Jk,k' , „• 1 

where Sj represents the impurity spin, S- = ^{Sx — i Sy) and is a 2 x 2 skew symmetric 
tensor. 

As before, we expand the NG spinors in their components centered around the nodes. 
Since we have four nodes the impurity Hamiltonian has terms which describe spin flip scatter- 
ing processes involving, in addition, eventual inter-node scattering processes. The strength 
of these scattering processes is determined by J{Q) where Q is the relative wavevector of a 
pair of nodes. There are four cases of interest: 

1. Q ~ 0, corresponding to scattering processes that do not mix nodes ("forward scat- 
tering"). The corresponding coupling constant is J(0) = Jq. 

2. Q ~ ^koii, which mixes nodes 1 with 3 and 2 with 4. This coupling constant is Ji. 

3. Q ~ 2A;oe2; which mixes nodes 1 with 4 and 2 with 3. This coupling constant is J2. 
For systems with exact tetragonal (square) symmetry Ji = J2. 
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4. Q ~ 2^0 (ci ± 62, which mixes nodes 1 with 2 and 3 with 4. These couphng constants 
are J^. For tetragonal systems they reduce to just one (diagonal) coupling J^. 

For example, consider an impurity seated at the Cu site at x = 0. As a crude approximation 
we may assume J{x) ~ J6{x). The Fourier transform tells us that all the couplings will be 
the same, and equal to J. A more realistic shape for J{x) would be a gaussian centered 
at the impurity site x = and decaying rapidly within a distance of the order of a lattice 
constant A. Thus we take J(x) ~ J/(27rA^) e^^^^'^^^''^^ will generate, for a generic fc-vector 
(which in our case will be 2A;oei, 2koe2, 2A;o(ei + 62) and 2/co(ei — 62)), 

J0) = J e--2^'''' (3.17) 

For the impurity at the origin in a tetragonal (square) lattice, all the coupling constants are 
real, with Jq > Ji = J2 > Jd- In the language of the fields introduced in Eq. ( |3.7| ) the 
impurity Hamiltonian now becomes 



Himp = |] J ^^^t .(p) J J^A,iiP 



a,b=l 



+ S- i: K+, i: J {^r,\. J ^<(-p') +h .c. (3.18) 



a,b=l i=l 



(27r)2^"'*^^^ ' J (27 

In Eq. (|3.18|) the indices a, b are the flavor indices which label the effective Dirac fermions 
species associated with each node. The indices i,j run through the spinor components (two 
per each NG spinor, i.e., per node) and label linear combinations of quasi-particles with spin 
up with holes with spin down. Also notice that p and p' now label small departures from 
the appropriate node. Using the fact that 2x2 is an antisymmetric matrix, we can rewrite 
Eqn. (|3.18| ) in the form 

u Q \- ^ r dp td9 t r 'dp' f^^dO' 

^ JO ZTT JO zvr JO ivr jo ztt 

a, 0=1 «=i 

^ + ^ dp f^^ d9 . , ^ ,dp' [^^ dO' t . ,x , . ^ 

a, 0=1 «=l 

The 4x4 matrices K^j, and K^^ used in Eqns. ( ^.18| ) and ( |3.19D are given by 



K 



ab 



Jo Jd Jl J2 ^ 

Jd Jo J2 Jl 

Jl J2 Jo Jd 

J2 Jl Jd Jo 



^ Jd Jo J2 Jl ^ 

^0 Jd Jl J2 

J2 Jl Jd Jo 

Jl J2 Jo Jd 



(3.20) 



The form of Eqn.( |3.lP| ) strongly suggests the following change of variables (particle- hole 
transformations) performed on the second component of all four flavors 

i^lM ^1,2(-P); ^2,2{P) V'2,2(-P); ^3,2(p) ^3,2(-p); ^4,2(p) ^4,2(-p) (3-21) 

to express the interaction term as a scalar product of two spin-| operators. 

We can now separate the modes and find an effective one-dimensional model. After 
integration over the angle variable 6, the fields involved in Eqn.( 3!T9D become 



dO 1 



(3.22) 



10 27T ^2 

Now we define, for each flavor a, an effective one-dimensional chiral (right moving) fermi 
field 



dla{p) 



d2a{p) = 



\p\ V'o,+,a(bl); for p > 0; 
bl ^o-,a(bl); for p < 0; 

bl V'-i,+,a(bl); for P > 0; 
/bl ^-i-,a(b|); for p<0; 

and Eqn. (|3l9|) can be recast as 

+ s- Y.Ktb r ^/h 4i(p) r ^/h 4(p') + h. 



(3.23) 



(3.24) 



'OO ZTT * J-oo ZTT * 



ab 



(3.25) 



Here we perform the change of variables suggested above by setting 



da2ip) ^ dl^i-p), for a=l,..,4 



(3.26) 



19 



and thus, in the definition given by Eqn. (|3.2^) we rename dia{p) as the d^a{p) component 
of an effective spin-i one-dimensional chiral fermion, and d^ai—p) as the dia{p) component. 
Please notice that this label is not equivalent to the spin of the original quasi-particles. In 
fact, the relation between these effective one-dimensional chiral fermions and the original 
quasi-particles is actually quite complicated. The (flavor) coupling matrices commute with 
each other (as required by the SU{2) spin rotation invariance) and can be diagonalized 
simultaneously by means of the following unitary transformation 



dL 



Uab du 



(3.27) 



where 2 =1 or |, the flavor indices a and b run from 1 to 4 and 



^1111^ 



u 



ab 



-1 1-11 
-1-1 11 
1-1-11 



(3.28) 



This rotation brings the coupling matrices Kj^j, and K^^^ to the diagonal form 



^ J[ ^ 









K 



^ J[ ^ 







ab 



-J^ 


-Ji 



040 



(3.29) 



with 



(3.30) 



J'l = Jo + Jd + Jl + J2 
J'2 = Jo ~ Jd + Jl — J2 
J3 = Jo + Jd — Jl — J2 
J4 = Jo ~ Jd — Jl + J2 

As one can see in Eqn. (|3.29| ) flavors 2 and 4 appear to have and Sy with the sign reversed. 
However this can be compensated by the following additional rotation in the spin components 
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(3.31) 



After all of these manipulations we find that the effective one-dimensional theory for this 
model is 



H, 



4 oo 
eff - 2^ 2^ ^ 

+ E E iJa/2) 
a=l a,u=^,l 



E{p)dl{p)dUp) 



Ztx * 



^ dp' f. — : , , , ' 

-^\^\p'\dav{p) 
oo ZTT * 



(3.32) 



In Eqn. (|3.32|) we dropped the primes in Eqn. (|3.27|) and in the effective coupling constants. 
The kinetic energy of the chiral fermions is -E(p) = \J vv'p. 

Eq. ( |3.32|) can be recognized to be exactly the non-marginal Kondo Hamiltonian that 



was discussed in ref. |29|. Hence, the effective Hamiltonian for a d-wave superconductor 
coupled to a magnetic impurity is essentially equivalent to a (multichannel) generalization 
of a non-marginal Kondo problem. There are four channels, one for each node. The channel 
degeneracy is generally lifted by the inter- node scattering. In fact, Eq. ( p. 301 ) shows that in 
the absence of inter- node scattering {i.e., Ji = J2 = Jd = 0) the four flavors couple to the 
impurity with exactly the same exchange interaction strength = Jq (a = 1, . . . , 4). For 
a strictly tetragonal system the couplings are ordered in the sequence J[ > J2 = J4 > J^. 
Intuitively one expects the channel with the largest coupling to dominate the low energy 
limit. In the extreme limit in which all inter- node and intra-node amplitudes are exactly 
equal one finds that channels 2, 3 and 4 decouple and that only the remaining channel 1 
couples to the impurity. Thus, in this limit, the physics of the system is that of a single 
channel non-marginal Kondo problem. 

Given that these two seemingly different systems are actually equivalent, most of the 
results found in ref, 



carry over to this problem almost without change but with a new 

we 



physical meaning and processes, in particular including pair breaking effects. In |29 



found that there is a critical value of the exchange coupling constant Jc, above which the 
impurity spin is screened. We also found there that the critical value Jc was of the same 
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order as the energy cutoff, which here is the superconducting gap Aq. The reason behind 
the existence of a finite Jc is that the effective interaction between the impurity and the 
normal excitations is momentum dependent and that it vanishes at small momenta (see 
Eq. |3.32| ). However, the same momentum dependence makes the effective coupling grow 
arbitrarily large at large momenta. This last behavior is unphysical and it results from the 
approximations, which are accurate at small momenta only. This observation motivates a 
simple redefinition of the model with a finite, momentum independent, coupling at momenta 
larger than a scale of the order of Aq/ {2'n\/vv'). 

In ref. p9| we showed that a momentum-dependent coupling is equivalent to a change in 
the density of states (DOS) for a theory with a momentum independent coupling constant. 
The model of ref. [|29|, and the model discussed above, have a DOS vanishing linearly with 
the energy. We consider a modified model with the DOS 

, ^ for |e| < Ao 
p(e) = <; 2™ ' ' - (3.33) 

for An < lei < D 



2nvv' ^0 

where Aq is the size of the superconductor gap away from the nodes. This change in the DOS 
is equivalent to a saturation of the coupling constant at the momentum scale {2tx\/vv'). 

In other words, we are assuming a linear dependence of the DOS with the energy around 
the gap nodes, up to the energy scale of the superconductor gap. For energies higher than the 
superconducting gap Aq, the normal quasiparticles are, for all practical purposes, identical 
to normal electrons. In a realistic cuprate superconductor, the band structure is actually 
rather complicated. Nevetheless, we can take into account the contribution of these states to 
the physics by considering a flat fermion band characteristic of a continuum spectrum from 
Aq up to a bandwidth which works as a high energy cutoff. As we will see below, the 
contribution of these states can almost always be ignored but they will enter in our results 
in two important places: (a) by shifting (downwards) the critical value of the coupling 
constant J^. and (b) in the scaling behavior for "half-filled" impurities. The shift in J^. is 
quantitatively important and it results in a downwards shift of Jc from the nominal value 
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of the superconducting gap Aq. Hence, we will assume that Jc < Aq. This happens if the 
scales of Aq and D are reasonably well separated. 



IV. LARGE Nc THEORY 

In the previous section we constructed a model for a magnetic impurity embedded in a 
d-wave superconductor and showed that it is equivalent to a special non-marginal Kondo 
problem. In this section we solve this model in the large N^. approximation, where N^. is the 
rank of the symmetry group of the impurity spin. In the physically relevant situation N^. = 2 
{i. e. , spin one-half). Clearly, in this situation A^^ is not large. Nevertheless we expect the 
large Nc theory to give a qualitatively correct description. We now proceed with a brief 
summary of the large- A^'c theory [|l^ as adapted [^] to the physical situation described by 



the Hamiltonians of the previous section. 

In order explore the physics of this system we extend the symmetry from SU (2) (spin) to 
SU{Nc) and look at it within the large- A^^^ approximation. Notice that, unlike the Coqblin- 
Schrieffer model, Nc is not related to a magnetic impurity in a higher spin representation. 
Similarly, the four flavors of fermions originate from the nodal structure of the superconduc- 
tor and are not related to an orbital-degeneracy as in the multichannel Kondo problem in 
metals. Thus the problem we want to study has Nc = 2 "colors" . The number of "flavors" 
is Nf = 1 if there is node mixing and Nf = 4 in the abscence of inter-node scattering. 
However, there is a subtlety in the treatment of the impurity once symmetry is extended 
from SU{2) to SU{Nc). For the group SU{2), the lowest representation for an impurity 
is 5 = 1/2 . For SU{Nc) many more representations are allowed. For example, the fun- 
damental representation, which has dimension Nc, is constructed by occupying an Nc-io\d 



degenerate multiplet with a single "slave" fermion For general N^ with the exception 
of Nc = 2, this representation is not self-conjugate or, in other terms, it is not particle-hole 
symmetric. Other representations can be constructed p6| by occupying the multiplet with 



Qf slave fermions. For Qf = Nc/2, which is available for Nc even, particle- hole symmetry is 
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exact. We will see below that particle-hole symmetry (self-conjugation) is a case of special 
interest. Notice that all choices of representation are, in principle, valid extensions from the 
physical SU (2)-invariant system. Similar caveats have to be made about the choice of a par- 
ticular generator in the algebra of SU{Nc) that will represent the Zeeman term for Nc > 2. 
In fact, in ref. [^] it was shown that some care has to be taken in this choice in order to 
describe a smooth weak-to-strong field crossover. In any event, we are only interested in the 
extrapolation of the results at A^c > 2 down to A^^c = 2 where there is no ambiguities but 
they are present for all Nc > 2. 

In reference p9|] it was shown that, after integrating out the fermion and impurity 
degrees of freedom, the impurity contribution to the effective action Sejf = /^-Flmp takes the 
form 



Firap = --J^Trln 

P (7=1 



Nf 



dr + ef + J2\(Pi\^Go{z) 



1=1 

'Nf 



' Nf 



dr\f\j:\<Pn-Qfef 



= Fin.p + / rfr I ^ I 5: |0P I - Q^ef I (4.1) 



^0 y/=i 

where 0; are the Hubbard-Stratonovich fields introduced to decouple the impurity in the 
large N formalism. The properties of the normal excitations is encoded in the function 
Gq{z) (where the complex number z = e + iX is the analytic extension of the energy). With 
the new definition of the DOS of Eq. |3.33| , the function Go{z), defined by 



oo 



G„(.) . - I -^f^ (4.2) 

Znv v' e — z 



now takes the form 



- In 



Ao 



for lei < Ao 

ReG'o(e + a) = ^ n u ^^3^ 

for lei > Ao 



-TT^ sgn(A) for lei < Aq 
ImGo(e + 2A) = <{ 2. ^ ^ ^ II ^^ ^^ 

- ^ sgn(A) for |e| > Ao 



We define 
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A 



TIV V 



(4.5) 



where Nj is the number of "flavors" . For the problem of the d-wave superconductor, Nj in 
principle is the number of nodes and Nf = 4. However, we showed above that the inter- 
node couplings are always different (and smaller) than the mtra-node coupling, this coupling 
anisotropy reduces to one the number of effective flavors. Hence, from now on, we will set 
Nf = l. 

At finite temperature T, the effective action of Eqn. ( ^?TD becomes an (infinite) series 
running over (imaginary) Matsubara frequencies. Using this approach the effective free 
energy becomes 



imp 



de 



2TTi J-oo 



In 



e^^ + 1 
— / de n{e) 6{e) 

TT J-oo 

where 6{e) is the phase shift ||3^, n{e) is the Fermi function 

1 



(e + a) + e; + (E,|0P)go(e + a) 
[e - a) + e; + (E, 101^)^0(6 -a) 



Explicitly we find 



tan 



n(e) 



A + 



e^^ + 1 



7r|e|A 
2 



eA In 



(A 0-) 



(4.6) 



(4.7) 



(4.8) 



Gq{z) has a branch cut and the jump of the function across this cut is energy dependent. 
This is an important difference with the usual Kondo effect in which the jump across the 



cut for the function Go{z) (see for example reference [|1^]) is energy independent and gives 
essentially the (constant) width of the resonance. This will not be the case any longer as the 
width of the resonance now becomes energy dependent. This marks an important departure 



from the "local Fermi liquid" (or the resonant level model) behavior [p!0| , ^ characteristic of 
the usual marginal Kondo systems. 

The large- A'^c analysis of this problem proceeds in the usual manner. Given the impurity 
free energy Fimp, a set of values of e/ and A that minimize this free energy are sought. The 
extremal values of e/ and A satisfy the saddle point equations (S. P. E. ) 
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BF- OF- 

^ = and ^ = (4.9) 
oA oef 

In the next subsection we will write explicit expressions for the SPE's and solve them. 

Thermodynamic magnitudes such as the impurity entropy 5'imp, the impurity contribution 

to the specific heat Cimp and the impurity contribution to the susceptibility Ximp as functions 

of temperature (and magnetic field) can be computed from the thermodynamic formulas 

OF- - d'^F d'^F 

^imp — Qrp ' "-^imp — Qrp2 ' -^»°1P ~ gjj2 K^'^'^J 

where, the total impurity contribution to the free energy Fjmp is given by 

v'iN 

-^imp — -^imp 

{A,ef,H,T) + TT ^ A - QfCf (4.11) 

Since both A and e/ are also functions of T and H, care must be taken to account for their 
contribution. However, since A and e/ satisfy the SPE's, we get 



dF OF- 

'^^ imp '^^ imp 



BF dF- 

J '-'^ imp '-'^ imp 



dT dT 

Thus, only the explicit dependence on T and H matters 



(4.12) 



V. SADDLE POINT EQUATIONS 

Using the formahsm of the previous section, the Saddle Point Equations (SPE) take the 
form 

1 85 
Qf = - den{e)—{e) (5.1) 
TT J-D oef 

and 

Jo J-D aA 

where D is the bandwidth cutoff. In general we will be interested in the regime T,H << 
Aq < D. In this regime, the contributions to the SPE from energies higher than Aq can be 
well approximated by setting T — H — 0. This amounts to setting the Fermi function to be 
n{e) ~ 0, for Aq < e < D, and n(e) ~ Nc, for —D < e < — Aq. The SPE's thus are a sum 
of two terms, one coming from energies |e| < Aq and one from I? > |e| > Aq. 
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A. Impurity Occupation 



The Saddle Point Equation Eq. ^]T] reduces to 

Qf = 



f|£|A 




(e + e A In 


e 





2 



1-^arctan ^ 



(5.3) 



For the reminder of this paper, we will be interested in the physics of this system close to 
the critical coupling constant. In that regime, the singlet amplitude A becomes very small 
and the asymptotic behavior of the SPE's in this domain can be evaluated explicitly. 

Thus, close enough to the phase transition, where A is very small, the contribution from 
the last term in |5.3| becomes 

/ Ao 



lim 



l-^arctan ^ 



A-A 



fAoA 



+ . . . 



(5.4) 



\Ao + e/y 3n \ Aq + ej ^ 
provided e/ + Aq > 0. 

For the remainder of this section we will consider the Saddle Point Equation in the case 
T = H = 0. In this case, the SPE takes the form 
^ A^^ r°° dz 



TT Je 



ttA 

2 



(5.5) 



ei/A z (^^^ +(Alnz + z/e-VA^)2 

where z = e^l^ Aa jy = jg clear that the integrand in Eqn.( |5.5| ) shows a crossover 

behavior at 



A In ^0 = 



-l/A 



Zq. 



In the regime v << A << 1, this equation has a root at large z, given by 

z, ^ e^/^) In e^/^) > e'/^ 

Taking into account the change in the behavior of its denominator, Eqn.(|5 
written in two pieces with the asymptotic form 

1 



2 



1 



+ 



Aln^ 



1 + Alnf (l + Alnf 
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+ . . . 



can be re- 



(5.6) 



and 



In 



1 + 



V 2 



(l + Alnf 



+ 



(5.7) 



Getting everything together, Eqn.(|5.1|) reduces to the expression 



Qf 



A A/-, A In 4 



+ 



21 + 1^ 21 + Aln^ 



+ 



Nr 



A\ni 



1 + A In ^ 



TT^A 



1 + 



+ ... 



(5.8) 



It is important to stress that, regardless of the approximations made in evaluating the 
integrals, the SPEl of Eq. |5.1| is a relation between A (the amplitude of the singlet) and 
the impurity Fermi level (in units of the gap Aq) ly at fixed occupation Qf. This relation is 
independent of the coupling constant and it must be solved first. For the problem that we 
are discussing here, the relation between u and A is singular, as implied by the logarithmic 
singularities in Eq. |5]^. We will see below that, due to the presence of this singularity, 
the impurity Fermi level e/ is no longer simply related to the singlet amplitude A. This 
phenomenon does not occur in the conventional Kondo problem in metals, where the DOS 
is constant. It occurs for systems with a DOS vanishing faster than linear with the energy. 
In this sense, the case of a linear DOS is a marginal system. 

Let X be the impurity filling fraction, x = Qf/N^. The solution of Eq. ^.8| takes the form 

1 2x 1 



z/(x. A) = \/e exp ( — 



+ 



[l + 0(A,AlnA)] 



(5.9) 



A l-2a; + A l-2a; + A. 
where e = 2.7172. . .. 

Hence, for generic values of x = Qf/Nc, the impurity Fermi level u depends on the 
singlet amplitude A through an essential singularity of the form exp(— const/A). As 
Qf ^ {x I), there is a crossover in the functional form of u which now behaves 
like exp(— const /A^), which vanishes much faster as A approaches zero. It is interesting to 
note that if the contributions from the states with energies between —D to — Aq had been 
neglected altogether, u would have vanished identically at Qf = for all finite values of 
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A. Since at Qf = the hamiltonian has an exact particle-hole symmetry, it may appear 
that V = e//Ao should have to vanish exactly at this point. In fact, it does not vanish due 
to states whose DOS violate the strict linear behavior of the DOS at low energies. 



B. Equation of State 



Let us consider now the second SPE, Eq.( [5.2| ). This equation relates A (the amplitude 
of the singlet) to the coupling constant (once the relation between e/ and A is known). We 
will regard this equation as an equation of state. 

AtT = and H = the SPE2, Eq.(|;2D can be written as 

TT^wl f^o TT e(e + ef) 

— ' de- 



Jo Jo 2 



^) +(e + e/ + eAln^ 



+ / de- Ao \ ^' (5.10) 



In the second integral of the r.h.s. of Eq.(5.10) it is useful to perform the change of variables 



u = e//e while, in order to treat the first integral or the r.h.s. of Eq.( ^.lO ) we use again 
z = e^/^(Ao/e). As above, u = e//Ao and un = tf/D. We can write 



1 _ TT p dz (l + ^ze 



90 2 L/A z^ ^ + Alnz)' 



+ ^ r ^ (5.11) 



'2, 

where we have defined the dimensionless coupling constant Qq by 

1 1 'k'^vI 



go Ao Jo 

The second integral on the r.h.s. of Eq. (|5.11|) can be shown to give the leading contri- 



butions 



TT du \ + u TT , D 1 TT / Ao \ , 9, 

In—- ^\{l-7a^) + ... (5.12) 



— / 9 ~ — m v \ i- 

2Juu u (^^1^% (1 + ^)2 2 Ao l + a2 2 V D 
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where = (7rA/2)^ << 1 in accordance to the hypothesis that A is small in the regimes 
in which we are interested. The first integral of the r.h.s. of Eq.( |5.11| ) is treated in the 
Appendix as an example of the approximations used. Retrieving here the results of Eg. ( |A8D 
we write Eg. ( |5.1lD in the form 

1 A 



90 2 VA„ 



TT 



2 l + (7rA/2)' 



TT 



A/2y 



+ o (a^) + 



Now we define the critical coupling constant as the limit for A — > of Eg.( p.l3[ ) 



1 

9c 



IT 



In 



D 



TT 

2 



For small A we obtain the scaling eguation 



(5.13) 



(5.14) 



11 . TT /TtA 

= 7rA + - — 



90 



4 



+ 
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VI. IMPURITY MAGNETIC SUSCEPTIBILITY 

In order to consider the effect of a magnetic field in our model it is necessary to proceed 
with some care. In principle we need to go back to the original model to look into the 
effects of a finite H. In terms of the Nambu-Gorkov spinors the "free hamiltonian" Hq now 
becomes 

Ho = J2 '^'(k) [(e{k) + /i) rg - A{k)n - h] ^k) (6.1) 

k 

Here the magnetic field H is multiplied by the 2x2 identity matrix. The conseguence of the 
introduction of a finite magnetic field is thus, the generation of a finite relative shift in the 
zero point for the energy, but not necessarily a finite density of guasi-particle states within 
the nodes of the gap (however, see below). The eigenfunctions remain unchanged but the 
eigenenergies are shifted by H 



E = - H ±^e'^{k) + A2(A;) (6.2) 
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Thus, after the expansion in small momentum around the nodes of the gap, the two- 
dimensional spinors will disperse with 

E ^ - H ± VW^p (6.3) 

where p has been defined before, in the model without field. 

It is not difficult to convince oneself, by going through the (several) transformations 
involved in the reduction to the effective one-dimensional model that, at the level of the 
one-dimensional hamiltonian, the magnetic field enters as a true magnetic field coupled 
now, to the one-dimensional chiral fermions. Thus, 

""^ dp 



Heff = T. E JZ^i'^P - HT,)dl{p)dM 



+ E E 



J-oo * J-oo Ztt * 



dp' 



(6.4) 



The change in the kinetic energy, depending on the spin polarization, changes the form of 
the function Go{u;, H) 

dp \p\ 



00 27r a; — \Jvv'p + Hr^ 

1 /-^ |e| de 1 ,^ , 1 /-^ |e| de 1 , ,^ ^, 
= wi-oo (a; + i/)-e 2^' + ^^) + W loo (a; - //) - e 2 " ^^"^^ 

where rs represents an SU (Nc) diagonal generator having r elements with eigenvalue and 

Nc — r elements with eigenvalue —1. In what follows we will take r = which respects 

the H — > —H symmetry of the SU{2) theory. For general r, a particle-hole transformation 

is not equivalent to if — >■ —H. Bout for r = iVc/2 these symmetry transformations are 

equivalent. In other terms, for general r this magnetic field breaks both the H ^ —H 

symmetry and particle-hole. Notice, however, that the representation of the impurity is 

determined solely by the charge Qf and it is unrelated to r. In the presence of the field, the 

impurity level has an effective fiUing factor 2Qf/Nc. We will see below that Qf — Nc/2 is 

a special case. For the physical case Nc — 2 there is only one possible representation (i e. 

spin S = 1/2) which corresponds to Qf = 1 = Nc/2. For general Nc these two situations do 

not necessarily coincide. 
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These changes will be reflected in the phase shift defined in section |^ In presence of 
a finite field the impurity free energy F^^p can be written as 

AT roo 

i^imp = de n{e) [6{e + H) + 6{e - H)] (6.6) 

ZTT J~oo 



where 5{e) is given by Eq.(|]H) and 72(e) is the Fermi function. Eg. ( |6.6| ) is manifestly invariant 



under the transformation H —>■ — H. 

The magnetization and the susceptibility are given respectively by 

dF- d'^F- 

which take the form 

In the limit T — 0, the function approaches a negative Dirac ^-function localized at 
e = 0. In this limit we find 

Mi^p(0,/f) = -^[5{H)-6i-H)] (6.9) 

ZTT 

Now we can use Eg. ( |4.8| ) to write an explicit expression for the magnetization 

M(0, H) = tan-i ( ] - ^ tan"^ ( ] (6.10) 

It is easy to see that in the limit H << ej(0), one has 

This expression shows that the impurity contribution to the magnetization vanishes as 
H^lnH with H ^ 0. As expected, the impurity magnetization vanishes as the field goes 
to zero thus showing that the magnetic impurity has been screenined. However, in a con- 
ventional marginal Kondo system, the magnetization vanishes lineraly with the field. Here 
instead we find a faster field dependence . 

It can be shown, using similar arguments, that a general expression for the impurity 
contribution to the magnetic susceptibility is given by 
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Ximp(T, H) 
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(6.12) 



At zero temperature the susceptibility becomes 

'85 



Ximp(0, H) 



2tt 



8e 



85 



H 



(6.13) 



Thus, we find that the susceptibihty at zero temperature and at low fields ( ^ e/(0)) is 



((^)' +{H + HA\nj^y + e}^ - Ae) (h + HA\nj^ 



(6.14) 



It should be noticed that in all of these expressions, the quantities and A are functions 
of the field H, with a limiting value e/(0) and A(0) for T = H = found in section 
0. The quantity A(0) should not be confused with Aq. Hereafter we will set A = A(0) 
and ef = ej(0). The magnetic susceptibility obtained from Eq. (|6.11| ) agrees with the limit 
ej » H oi Eq. ( |6.14|) and gives 

Ximp ~ H\n^ + h(i-^) + ... (6.15) 



In the opposite regime, H <^T <^ tf, the susceptibility is 

Ximp(T,0)^2iV,ln2 Tln(^; 



(6.16) 



To summarize, we find that in the low field limit the zero temperature magnetization vanishes 
like H"^ \a{AQ / H) . However, in contrast with the conventional "Fermi liquid" behavior of 
the Kondo effect in metals, the magnetic susceptibility is also found to vanish in the low field 
limit as ifln(A/if) and at zero temperature. Hence, in this non-marginal Kondo system, 
the magnetic impurity is overscreened even for a single channel of fermions. In the low 
temperature, zero field regime, the impurity susceptibility has a Tln(Ao/T) behavior which 
again shows that the impurity is overscreened. 
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VII. IMPURITY ENTROPY ANS SPECIFIC HEAT 



We can estimate the impurity contribution to the specific heat in the hmit T << H, 
in the screening regime. Using the SPE and some straightforward algebra, the impurity 
contribution to the entropy is 



imp 



OF 
df 



.A 



(7.1) 



Using the scahng e = xT this is 



'imp 



2tx J- 



[S{xT + H) + 5{xT - H)] 



(7.2) 



In the limit T « H , we may expand the phase shift around the point x = to get 

/OO 
dx x^ ^ (7.3) 

where Ximp(0,-ff) has been obtained in the previous section. The impurity contribution to 
the specific heat is then, equal to the impurity contribution to the entropy in this limit. Our 
result shows that the impurity entropy vanishes at T — >■ and finite field, for T << H. 
The general form of the specific heat is 



a 



imp 



27rT ^ ^ 9e 



85 



27r J-oo (e^ + l)2 

In the regime H <^T <^ ef we find 

A2 



8e 



Ac 



xT-H 



(7.4) 



Ci^p(0,T)^9C(3)iV,— T^ln(-f) 



(7.5) 



where ({3) is the Riemann (^-function at 3 and it is a number of the order of unity. 

Using the results of Eq. ( |6.16|) and Eq. (|7.5| ), we can compute the Wilson Ratio for the 
regime H <^T <^ ej and find 

a^^{H,T) 9C(3) 



Tximp{H,T) 2 In 2 



(7.6) 
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It is interesting that the ratio is still finite in spite of the fact the both the specific heat and 
the susceptibility behave very differently than in a Fermi liquid. 
In the high field limit T <^ H <^ Sf the impurity specific heat is 



TT 



2 



Q^^{H,T)^N,- ^-j TH\n{^) (7.7) 



which obeys the relation 



a^p(H,T)^-TximAO,H) (7.8) 



This result leads to a new Wilson Ratio 

Tx^ra^H, 0) 3 ^ ■ ^ 

which is essentially identical to the Wilson Ratio for the Kondo effect in Fermi hquids. 

Hence we found that in the strong coupling Kondo phase, the impurity specific heat at 
low temperature and low fields behaves like T^ln(Ao/T) and T i71n(Ao/i?) depending on 
whether H <^ T or T <^ H. Only for T <^ if we find the conventional linear T behavior 
of the (zero field) specific heat of the (marginal) Kondo effect in Fermi liquids. Notice 
however that the slope 7 of the specific heat in this regime is field dependent and behaves 
hke H\n{Ao/H). However, in spite of these differences, we found that the conventionally 
defined Wilson ratio is still finite but it is different in both regimes. 



VIII. DISCUSSION 

In this paper we constructed a model for the problem of a magnetic impurity in a d^^-yi 
superconductor. We solved this problem using the large- A^c approximation and found that 
there is a quantum phase transition from a phase in which the impurity is nearly free to a 
phase in which it is overscreened. We estimated the value of the critical coupling constant 
Jc- We found that Jc could be both smaller or larger than Aq, the gap of the d-wave 
superconductor, but it is certainly smaller than the bandwidth D of the electrons that 
participate in the superconductivity. 
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This result agrees with recent work by K. Ingersent [^] on a related system. Ingersent 
used a wilsonian numerical RG approach and found that the critical coupling runs off to 
the cutoff (strong coupling) unless either particle-hole symmetry is broken (for the band 
fermions) or additional high energy states with a flat DOS were added. In the problem of 
the d-wave superconductor the former possiblity is excluded by the superconductivity itself 
but the latter is required since such states are always there. In any event there is no reason 
to require that Jc should be smaller than Aq. In fact, even if > Aq, the Kondo scale 
Tk does not track Jc and it is almost always smaller than Aq (in fact, quite a bit smaller!). 
The value of the critical coupling constant is non-universal and it depends on details of the 
high energy physics of the system. Thus, our approximations have emphasized the role of 
the nodes and replaced the states above Aq by a "flat band". Clearly, the solution of the 
saddle point equations with the full band structure of the hamiltonian of Eq. ^]T] will yield 
a different (possibly smaller) value of Jc- The same caveats apply to the numerical RG 
calculation of K. Ingeresent, in which a speciflc discretization of the effective model is used. 
In fact, in most of his work, Ingeresent uses Wilson's logarithmic discretization which is very 
accurate for the Kondo problem in metals since it is tailored to reproduce the logarithmic 
singularities at high energies of the conventional (marginal) Kondo problem. In the case 
that we examine here, the system is very far away for its "lower critical dimension" . This 
approach should overestimate Jc, probably by quite a bit. In any event, the actual value 
of J itself depends on microscopic physics of the cuprates and there is no reason to believe 
that it should be tied to Aq. 

We investigated in detail the thermodynamic behavior (impurity susceptibility and spe- 
ciflc heat) in the overscreened phase where we found that, the impurity susceptibility van- 
ishes like H \n{Ao/H) (for T <^H <^Tk)oiT ln(Ao/T) (for H <^T <^Tk) with a crossover 
aX T H . For a Fermi liquid, the impurity susceptibiliy approaches aconstant value ast 
T ^ 0. The speciflc heat , on the other hand, was found to vanish like TH ln{Ao/H) (for 
T <^ H <^ Tk) or T^ln(Ao/T) (for H <^T <^ Tk) ■ In a Fermi liquid it vanishes linearly 
with T. The change in the power law behavior is an extension of the earlier work by Withoff 
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and Fradkin p6[. The additional logarithmic singularity is an indication that r = 1 is like 
an upper critical dimension for the Kondo problem The interesting quantum critical 
behavior, accessible for T,H ^ Tk was not discussed here and will be the subject of a 
separate publication [pO| ]. 

There are several important effects that we have not included here. One is the effect of 
random potential scattering which, naively may induce a non-zero DOS at the Fermi energy 
Ep = 0. Even if this effect is there, the effective DOS N{E) is very small. In |ll7| , |T8| it was 
shown that N{E) ^ exp(— const. /w) (where w is the witdth of the distribution) and that 
the elastic mean free path is exponentially long, i ~ exp (+ const. /t/;). S illCG Sit J Q we have 
a transition from a state with a divergent zero temperature susceptibilty (Curie-like) to an 
overscreened state with vanishing susceptibility, the rounding effects of a finite (but very 
small) DOS should be a very small correction if the material is clean. A more interesting, and 
perhaps more important, effect that was not included here is the presence of explicit pair- 
breaking by the vanishing of the amplitude of the d-wave order parameter at the impurity 
site. This effect should give rise to interesting Andreev like processes which may well alter 
the physics of this problem. We will discuss this problem elsewhere |^U[. Finally, corrections 
to the Nc —* oo limit remain to be estimated. 
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APPENDIX A: ESTIMATION OF INTEGRALS 



The integral of Eqn.( pTTO[) is a good representative for the approximations made in treat- 
ing the integrals of section 0. The expression given in Eq. (|5.12 ) is easy to obtain. As for 
the first integral of the r.h.s. of Eq.(5.11) we have 



Z2 



e'/^ / ^ = K + I (Al) 

2 Je^/^z^ (Af) +(i.ze-i/A + Alnz)' 



where the splitting of the integral corresponds to the + sign in the numerator. At the value 
zq = e^/^^ In-Zo there is a crossover in the behavior of the denominator of the integrand. 
For z < Zq the leading term is Alnz; for z > zq we can keep the term u z e^^l^. Also zq 
can be approximated by 



^0 ~ 1 — e 

I' J \v 



Then we can split the integrals as 

I < = - ^ - ^ 2 (A2) 

The change of variables t = In ^ makes the integration straightforward to give 



I< = ^ 
A2 



2 . \ / 2 



arctan — In 2:0 — arctan — - 



2 1 + Aln(f lii(fei/A)) 



On the other hand 



>_ TT 1 /-"^ rfz 1 

- 2 Z/2e-2/A Z ^2+(el/AEA^2 ^ ^ 



which can be integrated by partial fractions to give 



Similarly we have 



I " = ^ f 7 (AS) 

2 (l + Alnf^ 



^ " 2 ' A2 A/A ' (^2)^ + ^^^^ 
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and 



K 



1 



eVA 1 ^ 

< P ~ -z/ 

u Zq 2 



1 



2 i/2 j^^ z^ ^2 ^ (^^i/A V zo 2 (i + Alnf 



(AT) 



Hence, again K ^, as it was the case with I ^, can be neglected since its contribution is 
at least of the o(z/) ~ o(e~^/^) << A. The leading contribution comes from Eqn.(^6|) 
which can be recasted in terms of an exponential integral function and produces a constant 
contribution and a linear term in A. The leading contribution for small A give 
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TT 
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2A 



0{v) + 



l + M/2r (l + (7rA/2f 
With these results, the equation of state, keeping only the leading order contributions is 
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